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We study the dependence of the Kondo temperature on the gate voltage in a strongly blockaded 
quantum dot with a small single-particle level spacing. We show that the dependence cannot be 
fitted to that of the Anderson impurity model with the gate voltage-independent level width. The 
effect originates in high-order tunneling processes, which make a dominant contribution to the 
exchange amplitude when the gate voltage is tuned away from the middle of the Coulomb blockade 
valley. 
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I. INTRODUCTION 

In a typical transport experiment a nanostructure, 
such as semiconductor quantum dot, 1-3 is connected via 
tunneling junctions to two massive electrodes. Such de- 
vices often exhibit a logarithmic enhancement of the con- 
ductance G with lowering the temperature T, 



G(T) = Go + 
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K 
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with temperature- independent coefficients Go and Gk- 
This behavior is a manifestation of the well-known many- 
body phenomenon, the Kondo effect, resulting from the 
interaction of conduction electrons with an impurity pos- 
sessing additional degrees of freedom (see Ref. [4-6] for 
a review). The associated energy scale, the Kondo tem- 
perature Tk, characterizes also the dependencies of the 
differential conductance on the applied magnetic field 
and/or source-drain bias, which show a similar logarith- 
mic enhancement. 7 

The Kondo effect develops when a quantum dot has a 
non-zero spin in the ground state, which is guaranteed to 
happen when the dot has an odd number of electrons. 
The number of electrons N is controlled by the elec- 
trostatic potential V g on the capacitively coupled gate 
electrode. 1,4 ' 5 In the regime of a strong Coulomb block- 
ade, N is close to an integer at almost any V g except nar- 
row mixed-valence regions, where adding an electron to 
the dot is not associated with a large penalty in electro- 
static energy. Observable quantities, including Go,Gk, 
and Tk in Eq. (1), exhibit a quasi-periodic dependence 
on the dimensionless parameter Af = V g /8V g , where SV g 
is the distance between the mixed-valence regions. 4,5 In 
terms of Af ', these regions are narrow intervals of the 
width A <C 1 about half-integer values of Af. 

In this paper we study the dependence of the Kondo 
temperature of a quantum dot Tk on the gate voltage 
Af. Surprisingly, the dependence Tk{AT) has not receive 
much attention in the literature (see, however, Refs. [8] 
and [9]). In fact, experimental data are often 2,3 fitted to 
the expression 
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originally derived 10 for the single-level Anderson impu- 
rity model. In Eq. (2) Ec is the charging energy and A_\f 
is the distance in the dimensionless gate voltage to the 
charge degeneracy point, 



Ajs = l/2-\Af-Af \, 



(3) 



where Af is an odd integer. Eq. (2) is applicable for Ajy 
in the range 



A< Aaa < 1/2. 



(4) 



Since the Kondo effect is a crossover phenomenon 
rather than a phase transition, a precise definition of Tk 
is somewhat arbitrary. In particular, Eq. (2), as well 
as Eq. (13) below, is based on the perturbative renor- 
malization group 11 and defines Tk up to a gate voltage- 
independent numerical coefficient of the order of unity. 
The choice of the coefficient does not affect the valid- 
ity of Eq. (1), which is applicable in the weak coupling 
regime of the Kondo effect T » Tk- The value of the co- 
efficient, however, becomes important in the strong cou- 
pling regime T <C Tk- A survey of various definitions of 
Tk for the Kondo effect with spin— 1/2 on the dot 7 can 
be found in Refs. [3] and [12]. 

Apart from Eq, which sets the overall scale, the de- 
pendence Tk{N) as given by Eq. (2) is completely char- 
acterized by a single dimensionless parameter A. Despite 
its simplicity, Eq. (2) captures the most essential qual- 
itative feature of Tk{A/): The Kondo temperature has 
a minimum in the middle of the Coulomb blockade val- 
ley Af = Afo- At this point Tk is exponentially small, 
Tk/Ec cx exp(-l/A). 

Moreover, although Eq. (2) is inapplicable in the 
mixed- valence region Ajv < A [see Eq. (4)], it yields 
an estimate of the energy scale in this regime. Indeed, at 
Ajv ~ A the exponential factor in Eq. (2) is of order of 
unity, while the prefactor is of the order of A, resulting 
in Tk ~ EqA. For the Anderson model -EcA coincides 
with the tunneling-induced width of the single-particle 
energy level in the dot T [see Eq. (18) below]. Since in 
the mixed-valence regime the Coulomb blockade is par- 
tially "lifted", the level width T indeed represents the 
true scale characterizing the low-energy properties of the 
system. 
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Unlike in the Anderson model, the single-particle level 
spacing in a quantum dot S is much smaller than the 
charging energy Eg- 4 ' 5 It was shown in Ref. [8] that for 
S <C Ec the dependence Tk {Af) differs significantly from 
that prescribed by Eq. (2). However, this result was ob- 
tained in the limit when the contacts between the dot 
and the leads are almost open. 5,13 In this limit the dot 
is in the mixed-valence regime at all values of Af. In 
other words, the limit considered in Ref. [8] corresponds 
to A ~ 1, which is incompatible with the condition (4) 
of the validity of Eq. (2). In this paper we show that, 
contrary to the widespread belief, Eq. (2) is inapplicable 
to sufficiently large quantum dots in the weak tunneling 
regime as well, even though Eq. (4) is satisfied in a wide 
range of gate voltages. 

It should be noted that experiments [2,3] are not 
performed in the regime of open dot-lead contacts. 
Indeed, in the open contacts limit the temperature- 
independent elastic cotunneling 4,5,14 contribution to the 
conductance Go [see Eq. (1)] approaches G(0)/2; hence, 
the temperature-dependent Kondo contribution enters 
Eq. (1) with a vanishingly small coefficient Gk = G(0) — 
2G < G(0) (see Ref. [15]). Accordingly, in this limit the 
conductance is almost unaffected by the Kondo effect. 
On the contrary, in experiments aiming at indisputable 
realization of the Kondo effect, such as those described 
in Refs. [2] and [3], G(0) « G K > G . 

The rest of the paper is organized as follows: In Sec- 
tion II we describe the model of a quantum dot coupled 
by tunneling to the conducting leads. In Section III we 
derive the exchange amplitude of the effective low-energy 
Kondo model both in the limit 5 ^> Ec, corresponding 
to the Anderson model, and in the limit 5 <C E c , corre- 
sponding to a quantum dot. The results are discussed in 
Section IV. 



linearized near the Fermi level, which corresponds to a 
constant density of states v. 

The second term in Eq. (5) describes an isolated quan- 
tum dot. We consider the simplest model, 4,5 

H d = J2e n dld ns + E c (N-Af) 2 . (7) 

ns 

Here N = ^2 ns d) ns d ns is the total number of electrons 
in the dot and Eq is the charging energy. The single- 
particle energies e n are characterized by a finite level 
spacing 6 <C Ec- The Fermi level corresponds to eo = 0; 
this level is singly occupied when the number of electrons 
in the dot N = (N) is odd. 

Finally, the last term in Eq. (5) represents the tunnel- 
ing between the dot and the lead, 

Ht = t Q J2c{ s d ns +R.c. (8) 

nks 

A description of the dot-lead contact in terms of the tun- 
neling Hamiltonian Eq. (8) is possible 4,5 when the di- 
mensionless (i.e., in units of 2e 2 /h) conductance of the 
contact is small, 

g = AirTo/S < 1. (9) 

Accordingly, the tunncling-induccd width Tq = itvt\ of 
single-particle energy levels in the dot, the single-particle 
level spacing 5, and the charging energy Ec form a well- 
defined hierarchy, 

r < S < E c . (10) 
III. THE EFFECTIVE KONDO MODEL 



II. THE MODEL 

We consider a strongly asymmetric configuration, 
when the conductances of the dot-leads contacts are very 
different; this simplification does not affect the results. In 
this case the lead with a weaker coupling to the dot plays 
the part of a weakly coupled probe, Eq. (1) remains in- 
tact, 4 and for the evaluation of the Kondo temperature 
it is sufficient to consider coupling to a single lead, 

H = H C + H d + H t . (5) 

Here 

ffc = £&c]U. ( 6 ) 

ks 

describes electrons in the lead. For a lateral quantum dot 
system formed by electrostatic depletion of a 2D electron 
gas at the interface of a semiconductor heterostructure, 1 ' 2 
it is sufficient to take into account only a single propa- 
gating mode per dot-lead contact, 4,5,13 and can be 



When the gate voltage is tuned away from the mixed- 
valence regions, at Af close to an odd integer Ao, the dot 
has an odd number of electrons N w Ao and its ground 
state has spin S = 1/2. The low-energy excitations of 
the hamiltonian Eqs. (5)-(8) are then described by the 
effective Kondo model 4 

H = H C + Vp + J(s-S), (11) 

where p = J2 kk >s c L c k>s and s = T,kk>ss< c L(°W2)c fe , s , 
represent the local particle and spin densities of con- 
duction electrons, and the spin S is the projection of 
S = J2nss' d-tisi^ss' /2)d ns , onto the ground-state multi- 
plet of an isolated dot. 

The potential scattering term in Eq. (11) is responsi- 
ble 9 for the deviations of N from Ao, 

N-M,n-2vV. (12) 

The reduction of the original model (5)-(8) to the effec- 
tive Kondo model (11) is possible for | N— Af \ -C 1, which 
results in the restriction (3) on the allowed values of Af. 
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The exchange term in Eq. (11) leads to the Kondo 
effect characterized by the Kondo temperature 11 



T, 



K 



D (vJ) 1/2 exp(-l/z/J). 



(13) 



Here Do is the high-energy cutoff; it corresponds either 
to the threshold for the intradot excitations 5 or to the 
energy cost, 



E ± =2E c \Af-Af TlM, 



(14) 



for adding/removing an electron to/from the dot, 
whichever is smaller: 



D = mm{6,E ± } = min{<5, 2E c Ax). 



A. Anderson model (S ^> Ec) 



(15) 



We discuss first the limit <5 » Ec, corresponding to 
the Anderson impurity model. 10 Although, in view of 
Eq. (10), this limit does not correspond to a realistic sit- 
uation, it leads to the qualitatively correct dependence 
Tr-(A/"), see the discussion above. At S » Ec all but 
n = energy levels in the dot are either empty or doubly 
occupied. Projecting these levels out, we write 

N -> n^ + ni+Afo-1, 

where and % are the spin-up and spin-down occupa- 
tions of n — level. Substitution into Eq. (7) yields, up 
to a constant, 

H d = 2E C [n t rH - (Af - Af + l/2)(n t + n+)] . (16) 

The tunneling Eq. (8) induces transitions to states with 
n = level being empty and doubly-occupied. The tran- 
sitions are virtual and can be accounted for perturba- 
tively. In the second order in the tunneling amplitude t 
one finds 10 ' 16 Eq. (11), with 



J = 2ti 



E, 



1 



V 



2 \ E> 



E. 



> (17) 



and with E± given by Eq. (14). Reduction to the Kondo 
model is justified when \N — J\f \ <C 1. With the help 
of Eq. (12) and the second equation in (17), this condi- 
tion translates into Eq. (3) with the width of the mixed- 
valence region given by 



Ao 



4 Fp_ 

7T E C 



(18) 



Comparison with the first equation in (17) shows that for 
the Anderson model 



A = vJ Q , Jo - min{j(A0} = J(Af ). 



(19) 



Using Eqs. (3) and (18), the exchange amplitude (17) 
is written as 



vJ = 



A 



4A^(1 - A*) 



(20) 



Equations (13), (15), and (20) then yield Eq. (2) above. 
Accounting for higher order in to contributions results in 
a correction to J in Eq. (17). The correction is small, 
A J / Jo ~ ^Jo <C 1, 10,17 and its effect on the value of of 
the Kondo temperature Tr is negligible. 



B. Quantum dot (6 < E c ) 

The Anderson model result for J = min{ J(A/")} [see 
Eqs. (19)] remains intact even in the limit 5 <C Ec, when 
the model is no longer applicable. The leading correction 
now reads 17 AJ /Jo ~ r /<5 -C 1. Although the correc- 
tion is still small, it is by a factor E c /S 3> 1 larger than 
that in the Anderson model and results in an increase 17 
of Tk in the middle of the Coulomb blockade valley by 
the factor C with InC = (AJ/J )/vJ - E c /S > 1. 

The effect of the higher order in tunneling contribu- 
tions to the exchange amplitude turns out to be even 
more dramatic when N is tuned away from the middle 
of the Coulomb blockade valley Af = Nq. Indeed, it well 
known 18,19 that when Af is close to the mixed valence 
region, say, at Af w Afo + 1/2, transitions between the 
two almost degenerate charge states of the dot result in 
diverging logarithmic corrections to the tunneling ampli- 
tude. The origin of these corrections is again the Kondo 
effect, with the two charge states playing the part of the 
impurity spin. 13 ' 19 

Following Ref. [19], we project out virtual transitions 
to the state with Afo — 1 electrons in the dot, associated 
with the energy cost w 2Ec >• E + . This amounts to 
the introduction of a high-energy cutoff in Eqs. (5)- (8): 
\£k\, \e n \ < 2Ec- The projected Hamiltonian can be cast 
in the form of an anisotropic two-channel spin- 1/2 Kondo 
model with the physical spin s representing the channel 
index, 19 



H = E e m^lp^sm + hr z f z + ^ (r+f_ + r"f + ) . 

spa. 

(21) 

Here the "bare" (corresponding to the bandwidth D = 
2Ec) values of the coupling constants are I± — 2to, I z — 
0. In terms of |JJ.) and representing, respectively, 
charge states with Afo and Afo + 1 electrons in the dot, 
the pseudospin operators in Eq. (21) are given by 

f, = l{\m\-\m\), t + =ti = \m\- 

The operators ip in (21) are the relabeled operators c and 
d of Eqs. (6)-(8), 



''Ps^p.a—'fl d n 



Accordingly, the single-particle energies e p _ a = 
are characterized by the pseudospin-dependent density 



-p,a 



of states = 1/(5, v^, 
density is given by t - 



= v. Finally, the local pseudospin 

^2pp'aa' ^spa (°"aa' /2) ^P sp ' a ', 
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where components of the vector a are the Pauli matrices 
acting on the pseudospin degree of freedom. 
The scaling equations for the model (21) read 



dl z 



= 11(1- h), 



dl± 



d( ^ ZJ1 d( 
where C = \n(2E c /D), and 
~ 1 



h-\(^ + ll)\, (22) 



K + n) J z, t± = (*w) 1/2/ -l 

are dimensionless coupling constants. Since SU(2) sym- 
metry is broken, renormalization generates also correc- 
tions of the type Xw a V'lpaVVa 7 *- These terms lead 
to small pseudospin-dependent corrections to the density 
of states, 20 which we neglect. 
Equations (22) yield 



where 



and 



h{D) 



ll(D)=^ + ll{D), 



J x (0) = ^/^J^S 



7 2 ln(2£ c /.D), D^Ec^Td, 



(23) 



[HD/Tc 



Here 



T c = lEcie^ 121 



(24) 



(25) 



is the ener gy scale for the charge Kondo effect 19 and D <~ 
Ec\/To/5 is the value of the cutoff at which the two 
contributions to I± in Eq. (23) are of the same order of 
magnitude, I Z (D) ~ 7. 

Equations (23) and (24) are valid in the weak coupling 
regime I Z {D) -C 1 and as long as the bandwidth D ex- 
ceeds both the single-particle level spacing in the dot 5 
and the addition energy E + . (For I z (D) <C 1 the Knight 



shift-like renormalization E + 
glected.) At smaller D, 



(1 — I Z )E + can be ne- 



D<D, =max{5,E + (Af)}, 



(26) 



the Hamiltonian is given by Eq. (21) with rcnormalized 
coupling constants. At E + <C £L, it is equivalent to 
Eqs. (5)-(8) with the high-energy cutoff D*, and with 
the tunneling amplitude to replaced by the gate voltage- 
dependent t, which satisfies 



Avt 2 {D^) =6l£(D t 



(27) 



[We neglect a weak potential scattering arising from the 
z-component of the exchange in Eq. (21).] 

Further reduction of the bandwidth from down to 
Do [sec Eq. (15)] can be carried out without regard to 



the presence of multiple energy levels in the dot as these 
levels have been already accounted for in the renormaliza- 
tion of the tunneling amplitude. Projecting out the dot's 
excitations with the help of the Schricffcr-Wolf transfor- 
mation, 10 ' 16 we end up with the Kondo model Eq. (11) 
with the exchange and potential scattering amplitudes 



vJ = 



2vt 2 (D,) 
E+ ■ 



V = - 



2E, ' 



(28) 



cf. Eq. (17). Repeating the arguments that led to 
Eq. (18) above, we find the rcnormalized width of the 
mixed- valence region, 



A = 



4 r 

7T Eq ' 



TTUt 2 {5). 



(29) 



Further analysis depends on the parameters of the sys- 
tem. The richest behavior is realized in the limit 



(5/Ec) 2 « r /s « 1. 



(30) 



In quantum dots formed by electrostatic depletion of 2D 
electron gas, 1 the left-hand side of Eq. (30) is controlled 
by the size of the dot L, S/Ec oc 1/L for sufficiently 
large dots, 4 ' 5 while the right-hand side is proportional 
to the conductance of the dot-lead contact [see Eq. (9)]. 
Experimentally, these quantities are tuned independently 
of each other, 1 ' 4 and the inequalities (9) and (30) can be 
satisfied simultaneously. 

Equation (30) is equivalent to \n(8/Tc) <C y/5/To, 
which is compatible with the assumption that <~ S 
belongs to the weak coupling regime of the charge Kondo 
effect \n(8/Tc) 3> 1. The renormalized level width T in- 
troduced in Eq. (29) is then given by 



(31) 



4 \n 2 (6/T c ) 



and satisfies r <C T <C S. 

Close to, but still well outside the mixed-valence re- 
gion, the exchange constant (28) takes the form 



vJ = 



A <c Aw <c S/Ec, 



(32) 



similar to that in the Anderson model at A^f -C 1 [see 
Eq. (20)]. Note, however, that A in Eq. (32) is much 
larger than its Anderson model value A = A given in 
Eq. (18). 

Further away from the mixed- valence region Eqs. (23)- 
(28) give 

vJ= -^-\l+(Ec/5)A ]n 2 Aj, « « 1. 

4A./V L J 

(33) 

The gate voltage-dependent logarithmic correction in 
Eq. (33) remains small for all Af. 

Stretching Eq. (33) beyond its domain of applicability, 
and taking into account that at Ajv = 1/2 Eq. (33) rep- 
resents only a half of the exchange amplitude in the effec- 
tive Kondo model (at this point E + m E- and transitions 
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FIG. 1: Kondo temperature Tk as function of the distance 
Aaa to the closest charge degeneracy point. The dashed lines 
correspond to single-level Anderson impurity models with two 
different values of the tunneling-induced level width, F and 
r . The solid line represents a quantum dot with a small level 
spacing S <JC Ec- 

between charge states with Ao and A/o — 1 electrons in the 
dot make an identical contribution to J), we recover the 
result of Ref. [17], AJ /Jo ~ r /<5. This strongly sug- 
gests that the corrections considered in Ref. [17] and the 
ones studied in the present paper have a common origin 
and serves as an independent check that the procedure 
outlined above indeed captures the dominant contribu- 
tions to the exchange amplitude. 

IV. DISCUSSION 

The gate voltage-dependent corrections to the ex- 
change amplitude, and, therefore, to the Kondo tem- 
perature Tk, originate in the strong renormalization of 
the level width. Indeed, close to the middle of the 
Coulomb blockade valley Tk{M) is governed by the An- 
derson model expression Eq. (2) with A given by its 
bare single-level value A ~ T n /Ec [see Eqs. (18) and 
(33)]. However, the width of the mixed- valence region 
A <~ T /Ec 3> Ao depends on the much larger renormal- 
ized level width T > T , see Eqs. (29) and (31). 

Similar to the Anderson model, the width T represents 
the energy scale in the mixed-valence regime, while the 
corresponding value of A parametrizes the dependence 
T K (J\f) close to but still well outside the mixed- valence 



region [see Eq. (32)]. Unlike in the Anderson model, cf. 
Eq. (19), there is no simple relation between A and vJq = 
min{j/J(A/')} ~ A ; the latter determines the value of 

min{7V(A0}. 

These observations imply that the dependence Tk{M) 
cannot be fitted to the Anderson model's result Eq. (2) 
with the gate voltage- independent A. Instead, Tk{N) 
interpolates smoothly between the corresponding curves 
for two different Anderson models, as sketched in Fig. 1. 
This behavior appears to be consistent with the results 
of the experiments. 2 ^ 3 

In the above derivation we assumed that the tunneling 
amplitudes t n are identical for all energy levels in the dot, 
t n = to [ see Eq. (8)]. In large quantum dots with chaotic 
motion of electrons the amplitudes t n are not only dif- 
ferent, but random and statistically independent of each 
other. 4 ' 5,21 Accounting for these mesoscopic fluctuations 
does not affect our results qualitatively. The main differ- 
ence is that To and T are now random, although the rela- 
tion r > To still holds. Whereas the bare level width Tq 
has a broad Porter-Thomas distribution, 4 ' 5 ' 21 the renor- 
malized width T, being a sum of many statistically in- 
dependent contributions, is expected to have a narrow 
Gaussian distribution. 22 Accordingly, the exchange am- 
plitude J (and, therefore, Tk) become less random as the 
gate voltage approaches the mixed- valence regions. Note 
also that for t n ^ t- n the function Tk{N) is no longer 
symmetric about the middle of the Coulomb blockade 
valley, T K (N) + T K (2Af - M). 

Finally, as discussed above, Tk(N) in quantum dots 
with almost open contacts also differs 8 from that in the 
Anderson model. In this case the conductance of the 
dot-lead contact is large, 1 — g -C 1, and 

{vJ)- 1 ~ (E c /S)(l - g)sm 2 (nA M ). 

Description of the crossover between the result of Ref. [8] 
and our Eqs. (28)-(33), applicable for 5 < 1, requires 
an understanding of an intermediate regime between 
the strong and weak Coulomb blockade. Despite some 
progress in this direction, 23 a detailed theoretical descrip- 
tion of this regime remains an open problem. 

To conclude, in this paper we demonstrated that the 
gate voltage dependence of the Kondo temperature of 
a quantum dot with a small level spacing is drastically 
different from that in the Anderson impurity model. 
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